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The Rayleigh model of nonlinear Brownian motion is revisited in which the heavy particle of mass 
M interacts with ideal gas molecules of mass m <C M via instantaneous collisions. Using the van 
Kampen method of expansion of the master equation, non-linear corrections to the Fokker-Planck 
equation are obtained up to sixth order in the small parameter A = \Jm/M, improving earlier 
results. The role and origin of non-Gaussian statistics of the random force in the corresponding 
Langevin equation are also discussed. 
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I. INTRODUCTION 

There has been renewed discussion recently of the 
Rayleigh model of nonlinear Brownian motion, stimu- 
lated primarily by findings of new qualitative effects gov- 
erned by nonlinear stochastic processes such as rectifi- 
cation of fluctuations [1-4], stochastic resonance p|, etc. 
In the Rayleigh model the heavy Brownian particle of 
mass M is immersed in an ideal gas of molecules of mass 
m <C M and interacts with them through instantaneous 
elastic collisions. The gas is assumed to be so rarefied 
that collisions of its molecules with the particle do not 
affect the distribution of incident molecules, and also 
re-collisions can be neglected. In lowest order in the 
small parameter A = ^Jm/M, the velocity distribution 
function f(V,t) of the particle satisfies the second order 
Fokker-Planck equation which can be obtained by the 
conventional method, i.e. truncating the Kramers-Moyal 
expansion of the master equation to the first two terms 
and evaluating coefficients integrating the corresponding 
Langevin equation This procedure is consistent if 

the random force in the Langevin equation can be treated 
as white (with negligible correlation time) and Gaussian. 
However these two assumptions may be justified strictly 
speaking only in the ultimate limit A — * 0. For finite 
A, the Kramers-Moyal expansion has generally an infi- 
nite number of non-zero terms, and more careful analysis 
is needed to study effects of higher order in A pj. In 
this case, as was first shown by van Kampen, the dif- 
ferential equations for the distribution function do not 
have the form of the second order Fokker-Planck equa- 
tion but rather involve derivatives of order higher than 
two 0, . Although these higher order equations are lin- 
ear in the distribution function, they are sometimes re- 
ferred to as nonlinear Fokker-Planck equations because 
the corresponding Langevin equations involve nonlinear 
corrections of higher orders in A to the linear damping 
force. These corrections may lead to interesting physical 
consequences, such as additional terms in the fluctuation 
spectrum and modification of decay rate coefficients 0. 
Also, going beyond the lowest order is often necessary to 
account for delicate fluctuation-induced phenomena such 
as directional drift in the absence of systematic forces. 



If the temperatures on the left and right sides of the 
Rayleigh particle are different, the particle undergoes di- 
rectional movement even when the pressure on both sides 
is the same [2-4] . The systematic average velocity of the 
particle can be calculated as a perturbative solution of 
the Langevin equation with nonlinear corrections to the 
damping force |3(. Alternatively, one can use the corre- 
sponding nonlinear Fokker-Planck equation or an equiv- 
alent set of equations for the moments [4| . 

To derive a differential equation for the distribution 
function f(V, t) to desirable order in A, one has to ex- 
tract an explicit A-dependence in the Kramers-Moyal ex- 
pansion, transforming it into the expansion in powers of 
A, often referred to as the van Kampen expansion (VKE). 
The method essentially relies on an assumption of certain 
scaling properties of transition rates in the master equa- 
tion. For the Rayleigh model an explicit expression for 
the transition rates is available and all coefficients in the 
VKE can be found analytically. Note that another pop- 
ular exactly solvable model of Brownian motion, namely 
that of a particle coupled to a bath of harmonic oscilla- 
tors, is degenerate in the sense that all nonlinear correc- 
tions vanish identically. 

The Rayleigh model may be generalized in many ways 
to study effects of asymmetry of the surrounding bath [2- 
4], possibility of non-canonical distributions |1C( . finite- 
range interaction |ll| . strong non-equilibrium fluctua- 
tions etc. To analyze these generalizations it would 
be of help to compare their predictions with the results 
of the original model. In doing that we have found that 
the VKE for the Rayleigh model often appeared in the 
literature in forms which are either not quite accurate or 
incomplete, and lead to nonlinear Fokker-Planck equa- 
tions with serious defects. For instance, these equations 
may have a non-Gaussian stationary solution even when 
the Maxwellian distribution is an exact solution of the 
original master equation. Violation of the equilibrium 
condition can lead to large errors, such as an exponen- 
tially overestimated activation flux in the barrier cross- 
ing problem |13| . As was noted by van Kampen in 0, 
the Maxwellian distribution must satisfy each term of 
the VKE for the Rayleigh particle separately. While the 
equation of order A 4 derived in 8] has the Maxwellian 
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stationary solution, the equations obtained later by other 
authors do not always posses this virtue 0, hj. The 
extraction of explicit dependence on A for the Rayleigh 
model is slightly more complicated comparing with other 
applications of the van Kampen method and requires 
some care. One purpose of this paper is to give an ac- 
curate derivation of the equations for the velocity distri- 
bution of the Rayleigh particle up to order A 6 . These 
equations have the Maxwellian stationary solution and 
preserve positivity of the distribution function when devi- 
ations from equilibrium are small. Another purpose is to 
discuss some subtle points relevant to the derivation, such 



as the origin and consequences of non-Gaussian statis- 
tics of the random force in the corresponding Langevin 
equation, which has not been clearly articulated in the 
literature so far. Although the assumption of Gaussian 
random force was criticized in the literature , it is still 
widely applied to derive the second order Fokker-Planck 
equation, often without proper justification. In the last 
section we discuss why the random force is approximately 
Gaussian in lowest order in A, while corrections of higher 
orders may be essentially non-Gaussian. This fact has 
important implications which go far beyond the specific 
model considered here. 



II. KRAMERS-MOYAL AND VAN KAMPEN EXPANSIONS 



The conventional starting point in the derivation of the Fokker-Planck equation, as well as equations of higher 
order, is the master equation for the velocity distribution function f(V,t) 

= JdV'{f(V' ) t)W(V'^V)-f(V > t)W(V^V')}. (1) 

It is convenient to express the transition rates W(Vi — > V2) as a function of the initial state V\ and the transition 
length AV = V% — Vi, that is W(V\ — > V%) = W(Vi\AV). Then the master equation takes a more suggestive form 

df(V,t) _ 



dt 



d(AV){f(V- AV,t)W(V - AV\AV) - f(V,t)W(V\AV)} 
= fd(AV){v(V- AV,AV) - *(V,AV)j, (2) 
where ^(V, AV) = f(V, t)W(V\AV). Making in the expansion 

1 / P) \ 71 

^(V-AV,AV) = ^(V,AV) + J2^[-^Vg V ) nV,AV) (3) 

n— 1 ' ^ ' 

leads immediately to the Kramers-Moyal expansion, 
with coefficients a n given by 



a n (V) = J d(AV)(AV) n W(V\AV). (5) 

If the random force exerted on the particle by the bath molecules is white and Gaussian, one can show that only 
two first terms survive in the Kramers-Moyal expansion, which therefore turns into the second order Fokker-Planck 
equation. This can be proved using another equivalent representation for a n , 

a n (V) = Km±((V(t + T)-V(t)r), (6) 

which can be obtained from J5j writing the transition rates W in terms of the transition probability P(V, t\V' , t + t) 
as W(V -> V) = lim T ^o £P(V, t\V, t + r). Then © takes the form 

aJV) = lim - [ dV'CV'- V) n P(V,t\V ,t + t), 

r^O T J 



which is equivalent to JSJ. 
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The expression (0 for a n is often more useful than JSJ) since it does not involves transition rates W, which are 
usually unknown. On the other hand the average ((V(t + r) — V(t)) n ^, appearing in ©, can be readily found 
integrating the Langevin equation 

V = A(V) + F(t), (7) 
where A(V) is the damping force and F(t) is the random force with negligible correlation time, 

{F(t)F(Q))=TS(t). (8) 

Integrating J7J one gets 

V(t + t) - V[t) = f +T dt'{A(V(t')) + F(t% (9) 



which for small t can be consistently approximated as 

ft + T 

V(t + r)-V{t)^A(V(t))T + j dt'F(t'). (10) 

This approximation corresponds to a coarse-grained description with a time resolution To much shorter than the 
characteristic time for the relaxation of the particle's velocity ry and much longer than the correlation time for the 
random force tf, t~f <C tq <C ry. The limit r — >• in © means actually r — > To. Therefore, in © one can neglect 
time dependence of A(V(t)), but not of the random force which evolves significantly within the integration range 
t~t >t f . 

Using the approximation I|1U|) one can readily find ((V(t + r) — V(t)) n ) and then, from Eq. ©, the coefficients a n . 
For a white and Gaussian random force only two first terms survive in the Kramers-Moyal expansion, 

a x = A(V), a 2 = T, a n =0,n>2, (11) 

which therefore turns into the second order Fokker-Planck equation 

The result that in the Kramers-Moyal expansion only two first terms do not vanish relies entirely on the assumptions 
that the random force in the Langevin equation is Gaussian white noise. The above approach does not involve any 
parameters controlling applicability of this assumption, and the range of validity of the Fokker-Planck equation is 
difficult to analyze. This difficulty does not arise in situations when transition rates in the master equation are known 
explicitly In this case one can find a n directly from Eq. (|5j) without appealing to the Langevin equation and therefore 
without making any assumption about statistical properties of the random force. Using this method one finds that 
the condition a n — 0, n > 2 does not hold invariably, and the Kramers-Moyal expansion contains in general infinite 
number of terms. Nevertheless, it is still possible to get for f(V,t) a differential equation of finite order analyzing 
dependence of coefficients a n on the small parameter A and neglecting terms of higher order in A. The Kramers-Moyal 
expansion in the form Q is not appropriate for such a perturbation analysis because the dependence on A is implicit 
in Assuming that transition rates have certain scaling properties with respect to A, van Kampen modified the 
Kramers-Moyal expansion transforming it into the form of the expansion in powers of A. Keeping terms up to order 
A" in this expansion, one obtains for /(V, t) the differential equation of order n. 

III. VAN KAMPEN EXPANSION FOR THE RAYLEIGH PARTICLE 

For the one-dimensional Rayleigh model the transition rate has the form 

W(V\AV) = V 4 |AV| hi (v + \e- 2 AV^ , (13) 
where Jm{ v ) is the Maxwell distribution for the gas molecules, v is the number of particle per unit length, and 



M + m V f + A 2 



= VttW (14) 
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The transition rate has the scaling property 

e 2 W(V\AV) = 9(V\e~ 2 AV), (15) 

where 

Hv\z) = ^\z\f M (v + ^y (i6) 

To extract an explicit dependence of a n on A it is convenient to re-write Eq. (JSJ as follows 

a n (V) = J d(AV)(AV) n W(V\AV) 

= e 2n J d(e- 2 AV)(e- 2 AV) n (e 2 W(V\AV)y (17) 

or using i(T5|) . 

a n (V) = e 2n [ dtCHV\t). (18) 



Then the Kramers-Moyal expansion takes the form 

^-eM-^'w™™}- (19) 

71—1 X / 

where 

a n {V) = J dtc$<y\ti). (20) 

In (|19|l the dependence on A enters in two ways. First, e 2n as a function of A has the form 

e 2 " = A 2 >„(A), (21) 

where 

ma) = (rr^)" = 1 - - A2 + ^f^A 4 + . . . (22) 

Second, the expansion (|19f) involves dependence on the velocity the particle V which is obviously a function of A. 
If we restrict ourselves to small-range fluctuations about the equilibrium state, it may be reasonably expected from 
the cquipartition theorem that the ratio of the particle's typical velocity to that of surrounding molecules is of order 
y/rn/M. This leads to the second scaling assumption 

V = \x (23) 

where the scaled velocity x ~ A . Then the next step in the extraction of the explicit dependence on A is the expansion 
of a n (V) — a n (Xx) near V = 0, 

a„(Ax)=5>(f)^ (24) 
p=o P' 

where al P ' is p-th derivative of a n (V) at V = 0. Physically this expansion reflects the fact that the heavy particle 
is much slower than surrounding light molecules. The first term with p = corresponds to the particle which does 
not move at all (approximation of the infinitely heavy particle, A — > 0), while the next terms successively take into 
account the finite inertia of the particle. 

Substitution of the above expansions for e 2n and a n in (|19fl gives finally the desirable van Kampen expansion in 
powers of A 
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This form of the VKE is slightly different from that one usually finds in the literature and has the advantage that the 
dependence on A is entirely contained in the product A" +J V„(A), whereas the coefficients are A-independent. 
The VKE 1)25(1 can be written in the form 

k=l 

where the first six coefficients are 

Si(x,A) = -(A-A 3 + A 5 + ---) x (27) 
{a< 0) + Aa«x + ^A 2 a<V + ±\ 3 a[ 3) x 3 + ^A 4 a< 4) x 4 + ^ 5 af } x 5 + •••}, 

&(x, A) = ^(A 2 - 2A 4 + 3A 6 + • • • ) x 



{4 0) + Aa«x + -jjA'aJV + + + •■■}, 

S 3 (x, A) = -i(A 3 - 3A 5 + • • • ){4 0) + Aa^i + ±A a a< 2 > x 2 + |a 3 4 3 V + •••}, 

S 4 ( X) A) = i(A 4 - 4A 6 + • • • ){a 4 0) + A«5 + ^A 2 a 4 2) x 2 + •••}, 

S 5 (x, A) = -i(A 5 + • • • ){4 0) + Aa^x + •••}, 

S 6 (x, A) = I(A 6 + ---){4 0) + ■••}■ 

It is usually more convenient to write the result in the form of the expansion in powers of A, collecting in l|27[) terms 
of the same order, 

df(x,t) = f^ xkDkf{xt) (2g) 



dt 



k=i 



The first six differential operators are 
1 to' 



Dr = -^l, (29) 



n (i) a i (o) a 2 

ox 2 c^x 
n (o) d 1 (2 ) 5 2 1 (i) 9 2 1 (o) 9 3 

" 3 = ai to "2^ to" + 2 fl2 tox"""^ 3 av 

n (!) d 1 (3) 9 3 (o) d 2 1 (2) <9 2 2 1 (i) d 3 1 (o) a 4 
^ 4 = ai to""^ to" ~" 2 to+l" 2 to" "6" 3 to" + 24 a4 9x4' 

n (0) 9 1 (2 ) d 2 1 (4 ) 9 4 (i) <9 2 1 (3) C^ 2 3 1 (o) a 3 

to + 2^ to" "24 ai to" "" 2 to" + iT* 2 to" + 2" 3 to 



1 ( 2 )_^!_ 2^1 _J_ (0)J?_ 

12 ° 3 a^ 3 " + 24 a4 ax 4 " 120 as ax 5 ' 
n (i) 9 , 1 (3) a 3 i (5) a 5 3 (o) a 2 1 (2) a 2 2 1 (4) a 2 4 

to" +6^ to" -120 ai to" + 2" 2 to~2" 2 to" + 48" 2 to" 

,i (i) 93 _ 1 (3) 9 3 3 _ 1 (o) a 4 1 (2) a 4 2 _ 1 (i) a 5 1 (0) a 6 
+ 2 tt3 ax 3 " 36 " 3 ax 3 " e" 4 ax 4 + 48" 4 ax 4 " 120 " 5 ax 5 " + 72o ae ax 6 ' 

These formulas are valid not only for the original Rayleigh model, but also for asymmetric models [2-4] when properties 
of the bath on the left and on the right of the particle are different. Of course, the explicit form of the coefficients a[ I> ' 
are different for different models. In what follows we restrict ourselves to the symmetric problem when the transition 
rate is given by Ea. (|13fl . In this case, according to l|20|) and (|16fl . one gets 

an{V)= V - / dZe\Z\fM(V + Z/2) 
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and therefore 



«1 P) = 4 I dZC\Z\f$(Z/2)- (30) 



Here the Maxwellian distribution for the bath molecules is 

1/2 / 1 „\ m 



/«(.)-(£) exp -i^ , „_ "» (31) 



J/2' 

One can observe that offi = if n + p is odd. Then the equation (|28[1 contains only operators of even indices 
D 2 k, which involve an' only with even n + p. These coefficients can be expressed as linear combinations of the 



andf$(Z/2) = £f M (v)\ v ^ /2 . 
erve that 
olve a„ c 

integrals h = d^ 2k |£| /m(£/2). Using the first three of them 



j 64 3/2 1024 5/2 24576 7/2 

V 27T V 27T V 27T 



one can calculate all coefficients ol^ appearing in the equation of order A 6 , 

(1) V(J T -1/2 



„(3) - 



-Ii = —ao 

8 

8 32 



L5z/<7 3 5iA7 4 ,"'T 

^ J = h H J 2 h = aa i/2 

16 128 



„(°) 



f2l ^cr W i /o 

2 4 16 

4 J = Ji h H ^3 = 6acr 1/2 , 

2 4 1 8 64 d 

(!) V<J T 1fi -3/2 

to = — — Y 2 = — 16aer ' , 



4 3) 

„,(0) 



1 2 h = -48acr" 1/2 , 

8 32 1 

^7 2 = 32aa- 5 / 2 , 



CP1 i/o 
*F = -— 1 2 + — r/ 3 = 160aa- 3 / 2 , 
4 lb 

4 1} = /3 = -384a.-*/ 2 , 

4 0) = = 768aa- y / 2 . (33) 



In these formulas a = -7==. 

V Z7T 



IV. FOKKER-PLANCK AND HIGHER ORDER EQUATIONS 



The results of the previous section allow to write the van Kampen expansion for the Rayleigh model in the explicit 
form up to order A 6 . Truncation of (|28|l to terms of order A 2 leads to the second order Fokker-Planck equation 

= A 2 D 2 /0M), (34) 

where 
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and a = m/hsT. It was shown recently that the equation (|34|) . first derived by Rayleigh, can be recovered within 
a more general approach expressing coefficients in the A-expansion in terms of correlation functions for the random 
force and then taking the Markovian limit . 

The next approximation is the equation of order A 4 , 



df(x,t) 

at 



X 2 D 2 + A 4 D 4 



with 



/2. 



dx" 



6 dx 



r) 2 

2 -- 3/2 #T + 

OX'' 



cp- 

a - 1 / 2 — x 2 
dx 2 



+ 3 dx 3 + 3 9x 4 J- 

This equation is equivalent to that obtained by van Kampen in 
The equation of order A 6 reads as 



df(x,t) 



dt 



X 2 Do + X 4 D 4 



A 6 ZA 



}f(x,t), 



where 



Da = 



8;/ 



/2. 



dx" 
ft 2 
dx 1 



6 



1 



, 2 d_ 
dx 
d 2 



1 

120 



aV 2 £-x 5 



16 -5/2 d 4 

-Y a — 



8 

10 

T 



1/2^^4 _ 8(J -3/2 



fl 4 



dx 2 

5/2 j> 

&e 4 ' 3 " dx 4 
This equation appeared previously in an incomplete form in 



dx 

d 3 

dx 3 



Za- 3 ' 2 



d 2 



16 



-5/2 



9 5 



dx 2 
d 3 



15 



dx & y 



r 







missing the term a^J^x in _Pg 



(36) 



(37) 



(38) 



(39) 



One can show [l7| that the stationary solution of the 
master equation with transition rates in the form (|13fl is 
the Maxwellian distribution which for the scaled velocity 
x = X~ 1 V has the form f s (x) — C exp (— ^crx 2 ) . This 
distribution does not depend on A and therefore must 
satisfy each term in the expansion (|28|l separately, 



D n f,(x) = 0. 



(40) 



One can immediately check that the above expressions 
for Z?2, -D4, and Dq do satisfy this condition. 

Eqs. l|3f)jl or (|3"8"|l can be solved using an appropri- 
ate perturbation technique, however it is often easier to 
handle the equations for the moments (x n ) 0,01- 

It is known that approximations of the master equation 
by a differential equation involving derivatives of order 
higher than two may lead to solutions which are not pos- 
itive definite 0, . The reason why the Fokker-Planck 
equation (|34|) preserves positivity of / is because the right 
hand side of the equation, and therefore the time deriva- 
tive ft, are always positive at the points where f(x,t) 
as a function of x has minima. As a result, the minima 
become less deep with time, the initially positive solu- 
tion remains positive for all times. This is not true in 
general for Eqs. (|3l)|l and l|3^|l involving ^-derivatives 
of order higher than two. These derivatives can be of 
any sign at extreme points, which in principle may result 



in negative ft at minimum points. However, since the 
terms with higher order derivatives are of higher order 
in A, one may expect that for sufficiently smooth ini- 
tial distributions the sign of the right sides of Eqs. i|36|) 
and (|38|l is determined by the term with f xx . For ex- 
ample, the only terms in the right side of Ea. H36l) which 
can be negative at minimum points are those involving 
X 4 o~ 3 l 2 xf xxx and A 4 er~ 5 / 2 f xxxx . They are smaller than 
the term \ 2 o~ 3 l 2 f xx , which gives the positive contribu- 
tion, by the factors of order \ 2 (x/x c ) and X 2 (xth/x c ) 2 , 
respectively. Here x c is the characteristic length of the 
distribution f(x,t) for a given t, and Xtn is the scaled 
thermal velocity of the particle, x t n = X^ 1 ksT /M . 
Recall that the expansion method we used implies that 
the system is close to equilibrium and that x ~ xtn ~ A . 
Under this condition the term with f xx dominates and 
determines the sign of the right hand side of Eq. (|36|l at 
minimum points, which guarantees preservation of the 
positivity of the solution. 



The expansion method may be implemented also for 
non-equilibrium fluctuations [12j but in that case the pos- 
sibility to introduce a stochastic variable x, which would 
be of order one for all relevant times, is less obvious and 
has to be justified a posteriori. 
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V. DISCUSSION 

If the white random force F(t) in the Langevin equa- 
tion V(t) — A(V) + F(t) is also assumed to be a Gaus- 
sian process, then the conventional procedure outlined 
in Section 2 leads invariably to the second order Fokker- 
Planck equation i|12|l . no matter whether the damping 
force A(V) is linear or not. On the other hand, the 
van Kampen method, which does not require any as- 
sumptions about statistics of the random force, leads 
to the second order equation only to the lowest order 
in the expansion parameter A, when the damping force 
is linear A(V) = —7V. In higher orders in A, when 
the damping force A(V) involves nonlinear corrections 

A(V) = -71V-73V 3 (where 73/71 ~ A 2 ), the VKE 

leads to equations with V— derivatives of order higher 
than two. This means that the approximation of a Gaus- 
sian random force is legitimate to lowest order in A but 
not to higher orders. It is therefore inconsistent to take 
into account nonlinear corrections to the linear damping 
force A(V) — —7V, and at the same time to assume that 
the random force is Gaussian. Failure to appreciate this 
point may lead to wrong conclusions. For example, for 
the Rayleigh model the Maxwell distribution is a correct 
stationary solution to any order in A. On the other hand 
the assumption of a Gaussian random force would lead 
to the Fokker-Planck equation H12[l which for nonlinear 
A(V) has a non-Maxwellian stationary solution. 

The Gaussian property of the random force is usually 
expected to hold interpreting the force as a result of many 
uncorrelated collisions, and appealing to the the central 
limit theorem. However, one has to keep in mind that 
a decomposition of the total force exerted on the par- 
ticle into a regular damping and "random" parts is a 
purely mathematical procedure, which in general can be 
performed with an appropriate projection operator tech- 
nique. The random force F(t), obtained in this way gen- 
erally can not be interpreted as a superposition of many 
"physical" forces. Moreover, it is not even a dynami- 



cal variable: its evolution in time is not governed by the 
Newtonian propagator exp(iL) with the Liouville opera- 
tor L, but by a more complicated "projected" propaga- 
tor. One can get some insight in properties of the random 
force expanding it in powers of A which has the form [TT| 

F(t) =F {t)+\ [ dhSit-t^Foih) (41) 
Jo 

+A 2 f dt x I dt 2 S(t-t 1 )S(t 1 -t2)Fo(t 2 ) + ■ ■ ■ 
Jo Jo 

where S(t) is a non-Newtonian propagator the explicit 
form of which is not important for our purpose here. 
The first term Fo(t) in the expansion l|42[l is a dynam- 
ical variable and has a well defined physical meaning: 
it is a force exerted on the particle fixed in space (the 
limit of infinitely heavy particle). Conventional qualita- 
tive reasoning to justify the Gaussian property is quite 
applicable for this term. Moreover, for a large Brown- 
ian particle, interacting simultaneously with many bath 
molecules, the Gaussian property of Fo(t) can be proved 
analytically To lowest order in A the random force 
is just Fo(t) and therefore is Gaussian. On the other 
hand, there is no reason to expect that the same ar- 
gument should work for the "unphysical forces" repre- 
sented in l|42|) by terms of higher orders in A involving 
the non-Newtonian propagator. These terms in general 
are not Gaussian. In a future publication we shall explic- 
itly evaluate correlation functions (F(ti)F(t2) ■ ■ ■ F(tk)) 
for an exactly solvable model with parabolic interaction 
suggested in This would allow to construct the ex- 
pansion similar to the VKE but expressing coefficients 
in terms of correlation functions for the random force 
rather than transition rates W. In lowest order in A it 
was done in |ll| recovering in the Markovian limit the 
Fokker-Planck equation i|34|l • 

This work was supported by a grant from the NSERC. 
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